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Abstract 



We try to use scale- invariance and the 1/N expansion to construct a non-trivial 
4d O(N) scalar field model with controlled UV behavior and naturally light scalar 
excitations. The principle is to fix interactions at each order in 1/N by requiring the 
effective action for arbitrary background fields to be scale-invariant. We find a line of 
non-trivial UV fixed-points in the large- N limit, parameterized by a dimensionless 
coupling. Nether action nor measure is scale invariant, but the effective action is. 
Scale invariance makes it natural to set a mass deformation to zero. The model has 
phases where 0(N) invariance is unbroken or spontaneously broken. Masses of the 
lightest excitations above the unbroken vacuum are found. Slowly varying quantum 
fluctuations are incorporated at order 1/N. We find the 1/N correction to the 
potential, beta function of mass and anomalous dimensions of fields that preserve a 
line of fixed points for constant backgrounds. 
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1 Introduction 



We investigate the naturalness concept of 't Hooft [T] applied to 4d O(N) scalar fields. If there 
are scalar particles very light compared to the Planck mass, it must be due to a symmetry. We 
observe that one non-trivial scale-invariant RG fixed point in quantum scalar field theory would 
be enough to make small masses natural. For, setting masses to zero would buy us symmetry 
under scale transformations. We try to implement this idea by developing a method due to 
Rajeev [2] for constructing a fixed point in the large- N limit. 

1.1 Background and motivations 

Many discussions of 4d quantum field theory begin with massive A</> 4 theory. However, this 
model most likely does not have a non-trivial continuum limit [31 IU [6] . Our first motivation 
is to answer the question: 'Can one construct a non-trivial 4d scalar field theory?' 

Our second motivation comes from particle physics. The 2004 Nobel prize in physics for 
the discovery of asymptotic freedom in QCD has reminded us about the physical importance 
of quantum field theories with well-controlled ultraviolet behavior. Indeed, Yang-Mills theories, 
which have a gaussian high energy fixed point are at the heart of our best models for both 
the strong and weak interactions. In equilibrium statistical mechanics of 0(N) magnets, the 
gaussian fixed point controls high energy behavior while the lower energy dynamics is governed 
by a crossover to the non-trivial Wilson-Fisher fixed point [7]. 

However, the situation in 4d massive A</> 4 theory, which is the simplest and currently favored 
(but experimentally unconfirmed) model for W and Z mass generation, is less satisfactory for 
two reasons. First, A0 4 theory is based on the gaussian infrared fixed point, but does not flow 
to any UV fixed point to control the high energy behavior. Perturbatively, interactions become 
very strong at a finite energy ~ mexp [^ff - ] where m, A are the parameters of the model in the 
infrared (Landau pole). This is in contrast with asymptotically free theories or theories based 
on an interacting UV fixed point which might (at least in principle) be defined in the UV. Both 
numerical [3JH] and analytical [H [5] calculations suggest that in the absence of a UV cutoff, the 
theory is 'trivialQ. Unfortunately, the non-trivial Wilson-Fisher fixed point in quartic scalar 
field theory in 4 — e dimensions merges with the gaussian fixed point in 4d. As a practical 
matter, the lack of a UV fixed point in A0 4 theory does not prevent us from using it as an 
effective theory with a cutoff or as a perturbatively defined model like QED, over a range of 
relatively low energies. 

The second problem with Ad A</> 4 theory, is the naturalness problenil. In QED, a small 
electron mass (compared to mpi^^) is natural because setting it equal to zero gives QED an 
additional symmetry, chiral symmetry, which is not broken by quantum effects. On the other 
hand, in massive Xcfr theory, setting m = makes the classical theory free of any dimensional 
parameter. But scale- invariance is broken in the quantum theory due to the absence of any reg- 
ularization and renormalization scheme that preserves scale-invariance, there is a scale anomaly. 
In the absence of any symmetry to explain a small scalar mass in the quantum theory, natural- 
ness suggests that the lightest scalar excitation, the Higgs particle should have a mass of order 
the Planck mass. A very large Higgs mass, however, leads to other problems since (perturbative, 

The renormalized coupling constant vanishes identically and correlations satisfy Wick's formula. 
2 See appendix [X] for an explanation with examples. 
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essentially tree level) unitarity would then be violated [5|. The perturbative unitarity bound 
from a partial wave analysis of scalar exchange in W-boson scattering is estimated to be of 
the order of 1 TeV. Moreover, the likely triviality of the continuum theory implies a 'triviality 
bound' on the Higgs mass, which is also of the same order [9]. There does not seem to be 
any non-perturbative cure for these problems, they also arise on the lattice pJJ H] and in other 
analytical approaches [6] . 

Another issue is that the 1 -loop correction to the square of the bare Higgs mass is quadratic in 
the momentum cutoff and leads to the fine-tuning problem. If a large cutoff is to be maintained, 
either the effective Higgs mass is of the order of the cutoff or the bare mass must be fine-tuned 
to cancel the radiative correction (at each order). We already mentioned the difficulties with a 
very large Higgs mass. However, the fine-tuning problem needs to be treated with some care. In 
a renormalizable QFT (eg. Xcj) ), the use of a cutoff as a regulator is a matter of convenience. 
The quadratic divergence is absent in dimensional or £ -function regularization or Epstein-Glazer 
renormalization. In the end of a calculation, one sends all regulators to their limiting values 
before making physical conclusions. From this standpoint, a quadratically divergent self energy 
is not a deficiency of the model. 

On the other hand, one may be of the opinion that the QFTs of particle physics are effective 
field theories that come with a physical large momentum cutoff at either the Planck or other 
scale where new effects render the standard model inaccurate. This is similar to ones attitude in 
some (but not all) condensed matter physics contexts, where the crystal lattice is physical and 
not merely a convenient regulator. In this view, renormalizability is expendable. The cutoff is 
not to be sent to infinity, and the fine tuning problem mentioned above is indeed present. 

In QCD or other renormalizable models that are based on a UV fixed point (gaussian or 
not) to control high energy behavior, there is no need for the 'effective field theory-physical 
cutoff' view point. For, these models self-consistently predict low energy behavior irrespective 
of what the physics beyond the standard model may be or at what scale it may kick in. But 
in models such as , which (despite being renormalizable) seem to be non-trivially defined 
only in the presence of a cutoff, the latter viewpoint cannot be ignored. This leads us to an 
important distinction between the naturalness concept and fine tuning. The former requires any 
model to have an extra symmetry to explain a small parameter. But fine tuning, for its very 
definition, needs the model to have a physical cutoff. Fine tuning is potentially an issue for 
non-renormalizable models or models which need a cutoff in order to be non-trivial (eg. A</> 4 ). 

There is an empirical relation between naturalness and degree of fine tuning. Suppose we 
dogmatically insist on working with a cutoff even though our model may be renormalizable. 
Assume also that the model is natural, i.e. gains a symmetry if masses are set to zero. Then 
it is often the case that radiative corrections are 'protected by the symmetry': self energies are 
only logarithmically divergent rather than as a power of the cutoff. An example is QED. Due 
to chiral symmetry, the renormalization of the electron mass 2 is not quadratically divergent as 
one might naively expect, but only logarithmically divergent (the quadratic divergence persists 
in scalar QED, which does not have chiral symmetry). 

Our conclusions are (1) it is worth looking for a non-trivial scalar field model based on a UV 
fixed point; (2) it is worthwhile to look for a symmetry that would ensure naturally light scalar 
excitations and (3) If (1) and (2) are achieved, it is less important to worry about fine-tuning 
in the presence of a cutoff. 
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Despite our conceptual criticisms of Xcj) 4 theory, as long as a light Higgs is discovered, it will 
be possible to use the model to predict scattering amplitudes at relatively low energies. Thus it 
remains the default mechanism for giving masses to the weak gauge bosons. It may turn out to 
be an effective description of some more intricate framework devised by nature. 

Many alternatives have been proposed. Most studied is the use of supersymmetry to ensure 
light scalars [10] . The challenge here is to break SUSY without introducing new naturalness prob- 
lems. Technicolor models try to realize scalars as fermion composites [11] . Little Higgs models 
|12j attempt to use the Nambu-Goldstone mechanism in a novel way. In the Coleman- Weinberg 
mechanism |13] . classical conformal invariance of massless A</> 4 theory is used to explain small 
scalar masses. But this is not a naturalness explanation in the strict sense since conformal invari- 
ance is broken quantum mechanically. An idea to use scale invariance in a manner analogous to 
SUSY and soft SUSY breaking, has been proposed [14"] . Other approaches include higher deriva- 
tive models to expel one-particle scalar excitations as asymptotic states [15]; and the possibility 
of the gaussian fixed point being UV with respect to some non-polynomial potentials [16J. 

Richter [T7] has criticized naturalness. However, our definition of naturalness is not quite 
the same as the one he uses (in particular, a quadratic divergence in self-energy is not unnatural 
by our definition). Moreover, his criticisms seem to have more to do with the large number of 
parameters in the MSSM, than with naturalness as we understand it. We hope the examples 
in appendix |A] will help to improve the dismal score he gives the concept! Finally, one hopes 
that experiments at the LHC will add empirical discrimination to discussions on naturalness 
and electro-weak symmetry-breaking. 

1.2 Reasoning and summary 

Motivated by these considerations, we investigate whether there is any resolution to the difficul- 
ties of Xcj) theory that does not require adding new parameters or degrees of freedom beyond 
those of the standard modell. Can we construct an interacting scalar field theory in 4d, which 
has good UV behavior and supports naturally light scalar excitations? We argue that one pos- 
sible scenario is that such a model should be built around a nontrivial UV fixed point. Existing 
work does not indicate (indeed, almost rules out) a non-trivial fixed point in the neighborhood 
of the gaussian fixed poinlQ, so we will look farther afield. 

The two key concepts of this paper are (1) naturally light scalars from scaling symmetry 
in the quantum theory and (2) non-trivial quantum effective actions via cancelation between 
'action' and 'quantum fluctuations'. 

To find non-trivial fixed points far from the gaussian one, we give up thinking of a QFT as 
being defined by a classical action. The reason the classical action is a useful concept is that it 
provides a first approximation to the quantum effective action in the % — ► limit. Indeed, the 
contribution of the path integral measure is suppressed in this limit. By contrast, in the zeroth 
order of the approximation method we propose, both the 'action' and the quantum fluctuations 
from the 'measure' are comparable. In fact, both are strictly infinite prior to regularization 
and neither is scale invariant if regulated. However, their combined effect at the zeroth order 
of our scheme is to produce a finite 1 -parameter family of scale invariant (quantum) effective 

3 In this paper we work in an approximation where gauge and Yukawa couplings vanish. 

4 Halpern and Huang [TS] have argued that there may be potentials with respect to which the gaussian fixed 
point is UV. But this scenario is quite different from what we propose. 
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actions when regulators are removed. It is the latter that is physically observable and defines 
the theory. Such a scenario is not amenable to analysis via the loop expansion, traditional 
perturbation theory or a weak field expansion in powers of the scalar field. In particular, we 
give up perturbative renormalizability, but require it in a non-perturbative sense so as to ensure 
predictive power. We hope this is excusable since despite its great success elsewhere, insisting on 
perturbative renormalizability in scalar field theory leads to a model that is unnatural, lacking 
in good UV behavior and most likely trivial. 

To embark on the formidable task of constructing a non-trivial fixed point far from the 
gaussian one, it will help to have a small parameter to expand in. We investigate whether there 
is a 4d scale-invariant interacting O(N) scalar field theory in the 1/N expansion. The case 
of interest in particle physics is N = 4, since in the absence of gauge and Yukawa couplings, 
the symmetry group of the standard model is 0(4). Our approach is inspired by work on 
Yang-Mills theory, where we have learned that the theory has two different "classical" limits 
in which quantum fluctuations of some observables are small. Both the h ^ and N — > oo 
limits (holding the other parameter fixed) can be used as starting points for studying the full 
quantum theory. Might the large- N limit be of use in constructing a scale-invariant 4d scalar 
field theory? This possibility was pointed out by Rajeev [2]. 

If scale-invariance could be maintained even after including quantum effects, then such a 
scale-invariant model would define a fixed point. The way we achieve this, is to pick an action 
which is not scale- invariant, so as to cancel the 'scale anomaly' from quantum fluctuations. In the 
large- N limit, we actually find a line of fixed points^ parameterized by A. A is the dimensionless 
coupling constant of a A(/> 4 -type term which is marginally irrelevant when considered around 
the gaussian fixed-point but whose beta function vanishes in the large- N limit when considered 
around the non-trivial fixed point. In the spirit of Wilsonian renormalization, we should study 
the most relevant deformation of the line of fixed points. This corresponds to a mass term. 
Setting mass to zero should be natural, since we would gain scale-invariance by doing so. This 
line of scale-invariant theories would be UV fixed points with respect to the mass deformation 
and thus ensure good high energy behavior. For naturally light scalars via scaling symmetry, 
it is sufficient to have one fixed point. It may well be that when 1/N corrections are included, 
scale invariance can be maintained only for on^fl value of A = Ao ■ This would be acceptable, 
since m = 0, A = Ao would be natural due to scale invariance at that point. It would be a bonus 
if scale-invariance reduces the degree of divergence of self energy in the presence of a cutoff, just 
as chiral symmetry does in QED. But as a matter of principle, this is not necessary as long as the 
model is renormalizable and non-trivial when regulators are removed. Even more importantly, 
we would like locality, causality and unitarity. Locality is likely to be a subtle issue. We must 
also investigate the low energy behavior as we flow away from one of the UV fixed points along 
a mass deformation. 

Our model is to be constructed order by order in the 1/N expansion by requiring that 
the effective action for arbitrary backgrounds be scale- invariant. The physical output of this 
procedure is the effective action, not the classical action. Our starting point is a 4d N + 1 
component Euclidean scalar field model whose partition function can be written formally in 

5 A line of fixed points has appeared previously in 4d QFT. In Af — 4 super Yang-Mills theory, ffyM i s believed 
to parameterize a line of quantum fixed points. In a sense, things are simpler in M — 4 SYM since Qym = is 
the trivial fixed point and the line of fixed points may be studied perturbatively. 

6 There may be no non-trivial fixed point when 1/N corrections are incorporated, then our scenario would fail. 

7 It is convenient to work with N + 1 components cj>o, cf>i, • ■ • <t>N since we will integrate out all but <po . 
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terms of an 0(N + 1) invariant potential V 



Z = y [D0] exp - i y d 4 x|i|V0i| 2 + W(- 



(1) 



It is analyzed by introducing an additional field a via a Laplace transform so that the action 
is quadratic in fa. is the Laplace conjugat^l of the 0(N + 1) singlet rj = (section 

[2]). The reason to work with a instead of fa is that in the large- N limit holding % ^ fixed, 
a has small quantum fluctuations, while fa continue to have large fluctuations. The functional 
integral over all but one of the fa is performed, leaving only b = and a as dynamical fields. 
Keeping b facilitates discussion of 0(N + 1) symmetry breaking. 

The goal of a QFT is the determination of the (quantum) effective action. In section [3?2l we 
show that in the large- iV limit holding h fixed, there is a one parameter family of finite and 
scale-invariant effective actions when regulators are removed. For background fields B(x), E(x), 
it is obtained in an expansion around a constant background E G . In zeta-function regularization, 

r (B, E) = i J d 4 x (VB) 2 + E5 2 + AE 2 - ^ {(X - S )n(A)(E - E ) + 0(E - S ) 3 | . (2) 

Here A = — V 2 /E D and 11(A) is a specific function determined in appendix lC.3[ E Q is arbitrar){^|, 
it is not a parameter of the theory. It appears merely because we are studying the theory around 
a constant background field. Tq(B,T,) depends on the regularization scheme. For instance, to 
relate two schemes for constant backgrounds, A must be shifted by a finite additive constant 
which we have determined for zeta function, cutoff (sec. I3.1.ip and dimensional regularization 
(appendix IB. 1|) . The principle that the effective action be scale- free leaves the dimensionless 
coupling A undetermined, and it parameterizes a line of fixed points. In section 13.41 we show 
that Tq is free of scale anomalies. Effective actions in QFT usually involve proper vertices of all 
orders and are very complicated, ours is no exception. What is perhaps unusual is that unlike 
in the zeroth order of the loop expansion of X(j) theory, the effective action of our model in the 
zeroth order of the 1/N expansion already involves vertices of all orders. This merely reflects 
that in the large- N limit, we have already incorporated all quantum fluctuations of fa ■ ■ ■ <pN ■ 

The potential V(4> 2 /N) that leads to this effective action is not physically well-defined. The 
same is true about the quantum fluctuations coming from the path integral measure. Unlike the 
effective action, neither V nor the measure is finite when regulators are removed. Furthermore, 
neither the quantum fluctuations nor V is scale invariant. They depend on a scale M, which 
however mutually cancels to give a finite and scale free effective action. V(4> 2 /N) is not in any 
sense an approximation to the effective potential i.e., T for constant backgrounds. Nevertheless, 
V does appear in intermediate stages of calculations and many people want to know what it 

I .12 

is. Its 'finite part' in zeta function regularization in the large- N limit, grows as V(jy ) ~ 
iog(l|||VM 2 jV) ^ or l ar S e l^/N (appendix |D]). We have not yet determined its behavior for small 
\fa 2 /N, though there is no indication of singular behavior. V(rj) is most easily expressed in 
terms of the Laplace transformed potential W{a) . At N = oo the finite part of W(a) in zeta 
function regularization is (e is the base of natural logarithms) 



8 a is valued on a contour C from — ioo to ioo. Within our approximations, C must avoid R . 
9 When E is a non-positive real number there are further divergences, which we have not yet treated. 
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a Q is the constant background value of the field <r(x) appearing in the effective action. Though 
Tq involves arbitrarily high derivatives, W$ does not, indicating a form of locality. It is interest- 
ing to know whether this remains true in other regularization schemes and after including 1/N 
corrections. 

Next, we deform this line of fixed points by adding a mass term —m 2 a to W(a) , which 
explicitly breaks scale-invariance. In section H~T1 we determine constant extrema of the large- N 

2 

effective action and find a phase b = 0, a = where 0(N + 1) is unbroken and another phase 
b = ±m, a = 0, where 0(N + 1) symmetry is spontaneously broken to 0(N) . In section FOl we 
determine the masses of the lightest scalar excitations in the 0(N + 1) symmetric phase 

m , 4v / 3vrm , . , 

M t = 7 = and M„ = — ■ (4) 



We are not yet able to study the broken phase; it occurs at a = 0. Though W(a) is singular at 
<7 = 0, the effective action Y may be regular there (the effective potential is), we do not know. 

In section [5] we attempt to preserve scale-invariance after including quantum fluctuations 
in a and b at order 1/N. Some quantum effects were already present at N = oo, due to 
fluctuations of fa, •••,(/> jv which we integrated out. Determining the effective action at order 
1/N is complicated, so we make some further approximations which we hope to relax later. 
In sec. 15.1.11 we assume that quantum fluctuations are slowly varying and determine the 1/N 
correction to the fixed point potential, Wo + j?W\ that ensures the effective action F(B = 
0, S) = To + ]yTi is finite & scale-invariant for general background fields H(x). But there 
could be further divergences not canceled by since it is independent of B. In sec. 15.1.21 

we drop the assumption B = 0, but assume that background fields are constant and derive a 
Callan-Symanzik renormalization group equation (RGE) for the effective action in zeta function 
regularization. Imposing the condition that the j3 function of A vanishes and using the previous 
result for Wi(S) leads to a unique solution for the /3-function of the mass and anomalous 
dimensions of b and a fields. Within our (drastic) approximations, the model is renormalizable 
and the line of fixed points is maintained at order 1/N . It remains to analyze the phase in which 
0{N + 1) symmetry is broken, and also couple our model to fermions and gauge fields. Finally, 
we wonder whether there might be a dual/holographic/string description of our scale-invariant 
model by analogy with the AdS/CFT conjecture |18j. 



2 Lagrangian and change of field variables 

Consider an N + 1 component real scalar field (f>i,i = 0,2, ■ ■ ■ , N in 4 Euclidean space-time 
dimensions. We require the action to be globally 0(N + 1) invariant. The partition function is 




The factors of N have been chosen to facilitate a meaningful large- N limit, eft 2 is short for 
J2iLo 4>%4>i ■ ^(af) wm ^ e determined by requiring the quantum effective action to be scale- 
invariant. We will find the effective action in an expansion around N = oo holding h fixed. 
In this limit, rj(x) = (j) 2 /N > has small fluctuations and behaves classically but, fa continue 
to have large fluctuations. We will integrate out all the fa except fa and study the model in 
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an expansion around N = oo, where r\ and b = 4>o/\N have small fluctuations. We impose 
r\ = cj) 2 /N via a delta function 



Z = J [D<f>] J™ [Drj]e-iJ d^{\v^+NV( v )} -Q S (Nr,(x) - 2 (x)). 

Now insert the integral representation 

5(Nr) - cp 2 ) 



da 



,a{Nr,-4> 2 ) 



C 2iri 



(6) 



(7) 



C can be any contour from — ioo to ioo since the integrand is entire. Up to an overall constant 
that cancels from normalized correlations, © becomes 



Z = J[D(j>] J°°[Dr{] J [D<r]e-s / ' d ^U\^ 



2 +a4> 2 +NV(r))-Ncrr}} 



(8) 



a is a Lagrange multiplier enforcing the constraint r\ = (p 2 /N . Though C is not the real line, Z 
is real by construction. We separate b = (fro/y/N in anticipation of integrating out <f>\, . . . , cf)^ , 



Z = / [Db] / [Dcf>] / [D V ] / [Z) ff ]e 



E l J Ii{( v ^) 2 + <T ^}+ Ar ( Vfe ) 2 + Ar,j '' 2 + Ar ' l/ ( j ?)~ Ar<Tr ? 



(9) 



Henceforth does not include <j>Q . Now, we reverse the order of r] and a integrals and 

observe that the r\ integral is a Laplace transform at each space-time point x, 



\Drf\e- {N/2h) I ^[VM-H _ e ~(N/2h) J <PxW(o) _ 



(10) 



Then 



Z = / [Db] / / [Da]e 



-{1/2%) J d A x 



Y,^ 1 ^4>i) 2 +°4'l)+N(Vbf+Nab 2 +NW(cT) 



(11) 



Reversal of a and r/ integrals is allowed if W(a) is non-singular along contour C. Now we 
reverse the order of the 4> and a integrals. The 4> integral is a gaussian 



J [D(f>]e- {1/2h) £f=i / «^(-V a +<r)* 



det 



-V 2 + q 
2-kTi 



-N/2 



(12) 



This integral converges if — V 2 +cr has eigenvalues with positive real part. Since —V 2 is a positive 
operator, this is ensured if SRcr > (though the answer possesses an analytic continuation even 
to a with negative real part, as long as it stays off the negative real axis). Thus we get (up to 
an irrelevant overall constant) 



Z 



[Db] / [Da]e 



-NS{b,a) 



where 



S(b,a) 



2h 



h tr log(-V 2 + a) + J d 4 x|(V6) 2 + ab 2 + W[a) 



(13) 



(14) 
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Here, W(a) is obtained from V{rj) via a Laplace transform (|10|) . at each x 

e -{N/2h)W(a{x)) = f°° e -(N/2n)(V(rj(x))-a(x) V (x)) d ^ x ^ ^gj 



Conversely, ^(77) is obtained from W(c) by an inverse Laplace transform along contour C to 
the right of all singularities of W{a) 

^_ e -{N/2n)(W(a)+a V ) = e ~{N/2h) V(«) _ ^gN 
C 27TZ 

cr is a dynamical field, it carries space-time derivatives and complicated self interactions. We 
are interested in correlation functions of the a and b fields, which are defined as 

(a( Xl ) ■ ■ ■ a(x n )b{ yi ) ■ ■ ■ b{y m )) = ^J [Db] J[Da}e- NS{ - b ^a( Xl ) ■ ■ ■ a(x n )b( yi ) ■ ■ ■ b(y m ). (17) 

Positivity of r\ = (f) 2 /N implies that cr(x) is valued on a contour C from the south pole to north 
pole of the complex plane. On the other hand, the contour of integration for the b field is the 
real line, a has dimensions of mass 2 while b has dimensions of mass. 



3 Scale-invariance of the effective action at N = 00 



The interaction W(a) appearing in the action (I14p is expanded in inverse powers of N 
W{a) = W (a) + ^Wi(a) + ^W 2 {a) + ■■■. 



(18) 



We do not assume analyticity of W(a) at a = 0, so we do not expand it in powers of a. W(a) 
is to be determined by the principle that the theory be scale-invariant at each order in 1 /N . Of 
course, the action S(b, a) is also expanded in inverse powers of N 



S(b,a) 
where So 



So + jjSi + ^S 2 + 



h tr log[-V 2 + a] + J (fx\(Vbf + ab 2 + W (<r) j 



S x = ~ J d A xW x (a), S 2 = ^J d 4 xW 2 (a) etc. 



(19) 



Wi,2,3— must be chosen to cancel divergences and scale anomalies coming from fluctuations in b 
and a while W$ is chosen to cancel those from fluctuations in <j)\ - ■ ■ (j>jsr ■ The possible choice(s) 
of Wo,i,2... define the scale invariant fixed point(s) just as ||9^| 2 defines the trivial fixed point. 
W n are not counter terms in the perturbative sense. For a given choice of W n that does the 
job, there may be interactions with arbitrary coupling constants we can add to W n and preserve 
scale invariance. We want this to be a finite parameter family. W n are not restricted to be of 
any particular form. For locality, we wish to avoid arbitrarily high derivatives of a . 

N and Ti appear differently in S(b, cr) . N — > 00 is a 'classical limit' in which b, a have 
small fluctuations and is governed by the action So(b, cr). Ti. — ^ is also a classical limit, one 
in which the original fields <pi have small fluctuations. It is governed by the original action 
/ d 4 x[| V<j)\ 2 + NV((J) 2 /N)] . These two classical limits potentially capture different features of 
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the full quantum theory for a given W(a) . There is a priori no reason for the two limits h — ► 
and N — > oo to commute. 7i tr log[— V 2 + cr] is a quantum correction to the action in the 7l — ^ 
limit. But it is part of the 'classical' action in the N — ► oo limit. 

A theory is scale-invariant if its effective action T (Legendre transform of the generating 
series of connected correlations, which generates all 1PI or proper vertices, see appendix [B]), is 
scale-invariant. Such an effective action defines a fixed point of the renormalization group flow, 
r is obtained by averaging over fluctuations in b and a and is defined implicitly by 



e~ Nr{B ^ = J [D0\ J[D,] exp [ - n[s(B + /?, £ + q ) - /jg - ? g 



(20) 



B(x) and are arbitrary background fields while j3 and q ('varsigma') are the fluctuating 

fields, b = B + P, it = S + ^ . T(B, S) is calculated in a series in powers of 1/N holding h fixed 

t(b, s) = r (s, s) + IriCs, s) + ^2 r 2(^, s) + • • • • (21) 

To zeroth order in 1/iV, the effective action 

h tr log[-V 2 + S] + j d 4 x j(V£) 2 + SS 2 + W (S)| . (22) 



r (i?,s) = 5o(s,s) = i 



Trlog[— V 2 + is divergent. We must regulate the model so that it is finite. Then we 

must pick Wo(S) (which will depend on the regulator) in such a way that when the regulator 
is removed, Tq(B,T,) is not just finite but also scale- invariant. Similarly, W\ is determined by 
the principle that T\ be scale-invariant and so on. Actually, we will also have to allow for wave 
function renormalizations, but these appear only at order 1/A r , see sec. [5j 

3.1 Effective action for constant cr at iV = oo 

Tr log[— V 2 + appearing in the large- N effective action ro(.B, S) (i22"j) is most easily eval- 

uated for a constant background S(x) = S G . This is a physically reasonable first approximation 
if space-time inhomogeneities are small. S(x) takes values on the contour C from — ioo to 
ioo. We anticipate needing to pick the contour to avoid the negative real axis, so S G will be a 
complex number that lies off the negative real axis. 

3.1.1 Momentum cutoff regularization 

In momentum cutoff regularization, 

tr log[-V 2 + £J = / d 4 p d A q~5{p - q)~5(p - q) \og{p 2 + S a ) = 6(0) [ d 4 plog\p 2 + EJ 

J\p\<A J\P\<& 

j dpp 3 log[p 2 + S ] J d0 4 . (23) 



l 



(2k) 

Here the space-time volume is J d^x = O = (27r) 4 5(0) and / dU^ = 2ir 2 is the 'surface area' of 
a unit 3-sphere embedded in 4d Euclidean space. Thus, 



Ml C A 

h tr log[- V 2 + SJ = ^ J d P P 3 lo §(P 2 + S °) 
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64vr 2 

nn 

64vr 2 



2A 4 log(A 2 + £ ) - A 4 + 2A 2 E D - 2£ 2 log(A 2 + E ) + 2E 2 log £ 
2 A 4 log A 2 - A 4 + 4A 2 E D - 2E 2 log A 2 ( divergent terms) 



+2E log E G (non scale invariant finite term) 

— S 2 (scale invariant finite term) + terms that vanish as A —* oo 



.(24) 



We must pick Wo(£) such that the large- N effective action (|22p is both finite and scale-invariant 
when A — > oo. In sec J3.2l we do this for general T,,B. Here we get an idea of the answer by 
requiring that Tq(B, £) be scale-invariant for constant £ = S . To this end, we can pick 



-h 



64tt 2 



2A 4 logA 2 - A 4 + 4A 2 S C 



2S 2 logA 2 



32vr 



:S 2 logS . 



(25) 



This is the 'minimal subtraction' choice. Our principle that Tq(B, £ ) be finite and scale-free is 
ambiguous. We could add to this choice of Wo(E ), any scale-free finite term of the form A£ 2 , 
where A is a dimensionless coupling constant. Other terms such as m 2 E G or terms proportional 
to a higher power of E Q would involve dimensional coupling constants and would explicitly 
introduce a scale into the theory. The general choice leading to a scale- free Tq(B, E d ) is 



W (E o ,A) 



-h 



64vr 2 



2A 4 log A 2 - A 4 + 4A 2 E - 2E 2 log A 2 



32vr : 



-S 2 logE + AE 2 . (26) 



In the large- limit, we have a 1-parameter family of RG fixed points, parameterized by A. 
The addition of m 2 E D corresponds to a relevant perturbation of one of these fixed points. The 
addition of c^E™ for n > 2 corresponds to an irrelevant deformation, since coupling c n has a 
negative mass dimension. In the sequel we consider the mass deformed theory, where Wo(E G ) 
is a two-parameter (to, A) family 



-h 



64vr 2 



2 A 4 log A 2 - A 4 + 4A 2 E D - 2E 2 log A 2 



32vr 2 



E 2 logS + AS 2 -TO 2 E .(27) 



Wo(T, ) has a branch cut along the negative E G axis, consistent with our expectation that E(x) 
is valued on a contour that misses the negative real axis. The corresponding large- N effective 
action ()22[) for constant backgrounds B Q , E is 



To(B , £ ) 



n 

2 



m 2 E + A 



64vr 2 



E 2 + E OJ B s 



(28) 



Having found a line of fixed points, are they UV or IR? The answer can depend on the direction 
in which we flow from a fixed point. In this case, all the above fixed points are UV with respect 
to the mass deformation. The physical reason is that as we go to higher energies, the ratio of 
to to the energy scale will decrease and the RG flow will tend towards the fixed point. Recall 
that the gaussian fixed point (massless free scalar field theory) in 4d is UV with respect to 
mass deformations, but IR with respect to the quartic coupling A</> 4 . So far, in our model, 
the analogue of the quartic coupling, AS 2 is exactly marginal. Note that E has a canonical 
dimension of mass squared, while <j) h as dimensions of mass. 

The presence of a 1-parameter family of UV fixed points is a godsend. It means that for any 
values of A and h, we can set to = and get an extra symmetry, scale-invariance. Thus, the 
mass parameter can be naturally small (at least at N = oo). An interesting question is whether 
the line of fixed points can be maintained after including effects of quantum fluctuations in a 
and b in a 1/N expansion (see sec. [5]). 
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3.1.2 Zeta function regularization 



We repeat the evaluation of the large- N effective potential by regularizing tr log[— V 2 + X ] via 
zeta- function regularization. This method directly prescribes a finite part for tr log[— V 2 + S G ] 
in the unregulated limit. The finite part is not scale- invariant. We will pick Wo(£ D ) to cancel this 
non-scale-invariant quantity, so that Tq(B, S ) is both scale-invariant and finite. This procedure 
means we do not need to prescribe how Wq(£ ) must depend on the regulator, but just the finite 
part of its limiting unregulated value. Such a short-cut is not possible in other regularization 
schemes such as momentum cutoff (sec. 13. 1 . 1|) or dimensional regularization (appendix IB.ip . For 
this and other reasons, we will use the zeta-function regularization in the rest of the paper. On 
the other hand, comparison of different schemes allows us to better understand what is scheme 
independent. We find that the presence of a 1 -parameter family of fixed points is scheme 
independent. Moreover, the finite and scale-invariant effective potential obtained in the three 
schemes are the same up to a finite relabeling of A . To define the finite part of tr log[— V 2 + E ] 
by zeta-function regularization, let 



C(s) = tr [-V 2 + SJ 



d 4 p d A q 5 A (p — q) 



\p 2 + 



n 



i 



d 4 P 

Wf b 2 +s r 



(29) 



Q = J d 4 x = (27r) 4 5 4 (0), appears because S(x) is homogeneous and we are taking the trace 
of an operator that is diagonal in momentum space. C( s ) is analytic for $ts > 2. By analytic 
continuation, we get a meromorphic function C(s) . Away from singularities, 



C'(s) = 

So, if C(s) is analytic at s = 0, C'(0) 
COO 1 f f°° p 3 dp 



tr 



(2tt)< 



log[-V 2 + S ] 
[-V 2 + S ]* 

tr log[— V 2 + S Q ] . Now let us calculate 

1 . o X 2 -* 1 



r2vr 



(p 2 + T, ) s (2?r) 4 2(s-1)(s-2) 16vr 2 (s - l)(s - 2) 



(30) 



(31) 



C(s) has simple poles at s = 1, 2 but is analytic at s = so, 

C'(0)_3Sg S 2 logE Q £ 2 log[£ e- 3 / 2 ] 
Q 64vr 2 32?r 2 _ 32vr 2 

log[e- 3 / 2 S ] 



64vr 2 32?r 2 
tr log[-V 2 + SJ 



E|0 
32tt 2 



(32) 



(j32[) is not scale-free due to the logarithm. The effective action at N = oo for constant back- 
ground fields B, £ is 



r (-B o , So) = — 



01og[e- 3 / 2 ^] + S o i? 2 + ^o(5] o ) 



(33) 



The parameter M with dimensions of mass sets the scale for the logarithm and breaks scale- 
invariance. The choice of Wq{g) that ensures Tq(B,T, ) is scale-free (for m = 0) is 



W (a o ) 



-m 2 a + Act 2 



32vr 2 g[ M 2J - 



(34) 



We added a mass term, a relevant perturbation away from the line of fixed points parameterized 
by A. We will prove in sec. 13.41 that the scale anomaly in ro(-B , S G ) vanishes for this choice of 
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Wo(a ) . We avoid cubic and higher powers of a Q as before. Note that the terms in Wq(u ) are 
of different orders in Ti but of the same order in jj . For this choice of Wo , we get 

T (B o , S ) = ^[-m 2 S + A£ 2 + E B 2 }. (35) 

M cancels out from the effective potential, which is scale- free for m = 0. Though Wo(a ) 
has a branch cut along the negative E Q axis, Fo(-B ,£ ) (and So(b ,a )) is entire. It is very 
interesting to know whether this is true for general backgrounds. The method of expanding 
in inverse powers of E Q that we use in the rest of this paper prevents us from answering this 
question here. Comparing with section 13.1.11 we see that irrespective of regularization scheme, 
there is a one parameter family of fixed points parameterized by A. However, the definition of 
the coupling depends on regularization scheme, 

-^zeta function -^cutoff (3471-2 ' (*^) 



3.2 N = oo effective action expanded around a constant background 



In section [3TT1 we calculated the N = oo effective action (|22[) for constant B,T,. Allowing for 
arbitrary backgrounds B(x) is easy, the difficulties lie in non-constant E(a;). Here, we get an 
expansion for Fo(-B, E) in powers and. derivatives of XI — S wher6 S is a constant background 
whose value is not a negative real number or zero. It would be interesting to calculate Tq by 
complementary methods too. From appendix [Cl in zeta function regularization, 



tr log[- V 2 + £(*)] = ^ log[£ e- 3 / 2 ] + J 



d A x 



32vr 2 bL ° 1 J 16tt 2 _ 

(s - s )n(A)(s - s„) + e>(£ - s ) 3 



S D log[E /e](E - E c ) + i(E - E ) 2 log E 



(37) 



This reduces to (j32|) for constant backgrounds (E = E Q ). The 0(E— E G ) 3 and higher order terms 
can be obtained by the method of appendixO But these higher order terms are scale-invariant. 
We will show in section [331 that the part of tr log[— V 2 + E(x)] that is not scale-invariant, is 
restricted to the first three terms on the rhs of (137]) . Note that A = — V 2 /E G and 



n(A) = £ 



(-Ay 



A(3A + 2) -2(A + l) 2 log(l + A) 



i n(n + l)(n + 2) 
Thus, the effective action at iV = cxd is 



4A 2 



A A 2 

"6 + 24 



60 



+ •••(38) 



r (B,S) = \Jd A x (VS) 2 + cTS 2 + ^ (S) + I ^2|iE 2 log[E o e- 3 / 2 ]+E o log[E o / e ](E-E o ) 



+\ log[E ](S - S ) 2 - (E - E )n(A)(S - S ) + £>(£ - E f 



(39) 



In deriving (|39|) we did not assume E is slowly varying but rather that E — E D is small and that 
E is not a negative real number or zero. 
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3.3 Fixing the interaction at N = oo by requiring scale-invariance 



Wo (a) must be chosen so that Tq(B, EJ) in (j39j) is scale-free. The choice that does the job is 

Wo(a) = Xa 2 - m 2 a - ^{ f log[^] + * log[^](<r -a ) + \ log[-^](* - a ) 2 } (40) 

with m = 0. For m / we have a mass deformation. M sets the scale for logarithms, but 
cancels out in the N = oo effective action 



/ ^ (V5) 2 + Si? 2 - m 2 S + AS 2 - |(s - S )n(A)(E - E c ) + 0(S - S ) 3 | .(II) 



ro(-B,S) is now free of divergences (for S Q ^ R~ ) and defines the large- N effective action for 
background fields S whose deviation from a constant S D is small. The higher order terms in 
<j = S — S a are all finite, scale-free and calculable by the method of appendix [Cj The fact that 
Tq(B, T,) is not quadratic when regulators are removed, indicates that our theory is not trivial. 



3.4 Cancelation of scale anomaly 



We show that the effective action Tq(B, S) in (|4ip is scale invariant for m = (m ^ is treated 
subsequently). Coordinates and fields are rescaled according to their canonical dimensions. M 
sets the scale for logarithms and explicitly introduces a scale into the theory. To encode this 
physical fact, M is not rescaled, the same is true of m. We define dilations as 



D a x^ 



~ l x^, D a b = ab, D a a = a 2 a, D a X = A, D a V = aV. 



(42) 



The generator of infinitesimal dilations 8q is defined as 

r D l+£ f-f 

un I — lim 

5dx^ = -a 
and may be represented as 



S D f 

5ob = b, 5do~ = 2a, 



6 D \ = 0, 5° D V = V, 5 D dx 



-dx 



Sd 



dx^ 



+ b(x) 



d 



db(x] 



+ 2a(x) 







da(x) 



(43) 



(44) 



We now show that Tq(B,H) (I22p is invariant under dilations for (not necessarily constant) 
backgrounds B(x) and S(x) if Wo(£) is chosen according to (f4T)j) with m = 0. First, 

D a J d 4 x{{VB) 2 + TiB 2 } = J d 4 x{{VB) 2 + £5 2 } 5 D J d 4 x{(VB) 2 + £5 2 } = 0. (45) 

VFo(5]) and tr log[— V 2 + E] are the only terms in (|22p with non-trivial (in fact inhomogeneous) 
scale transformations. Wq transforms as 



D a J d 4 xW (Z) 



d 4 xW (Z) 
d 4 xW (Z) 



d 4 x 



16vr 2 
hT? o 0, log a 
8^2 



^S 2 log[a 2 ] + S log[a 2 ] ? + \ log[a 2 ] ? 2 

1 + M + ifl 1 

2 + E„ 2S2 



14 



5 D J d 4 xW (Z) 



8ir 2 ■ - 



2 E 2E 2 , 



(46) 



where ? = E — E D , E D is a constant background and (/) = jj / cftxf. On the other hand, 

tr log[-V 2 + E] = -C'(0) and £(a) = tr [~V 2 + E]~ s 
=► D a C(a) = a- 2s C(s) and D C'(s) = -2C(s)a" 2s log a + a- 2s ('{s). 



(47) 



Now set s = and use the result for £(0) calculated in Appendix IC .41 



£> o C , (0)=C'(0)-2C(0)loga ^C'(O) = -2C(0) 



8vr 2 



1 , (?) , <? 2 > 



fofc tr log[-V 2 + E] 



-^dC'(O) 



~8^" 



2 + E D + 2E 2 , 

2\n 



i , (?) , (O 



+ + 



2 E 2E 2 , 



(48) 



We see that the scale anomaly of tr log[— V 2 + E] exactly cancels that of Jd xWq(S). Therefore 
the large- N effective action ([41 p (for m = 0) is invariant under scale transformations f|42j) : 
5£)Fo(-B, E) = 0. However, Tq with a mass term (I4ip is not scale- invariant. In fact, 



8 D T (B,-£) = -2 J (fx m 2 E(x) = J d 4 x m 2 E(x) 



(49) 



Define (3^ = m. Though Tq is not scale- invariant, it satisfies a 'renormalization group equation' 



(fo + /3S 1 ^)r (fl,E) = o. 



(50) 



ft? is the N = oo 'beta function' of mass. In general, j3 m = ffi 1 + ^/3f + + • • ■. We 

could also have /? A at higher orders in 1/N , but we want to ensure (3 X = for at least one value 
of A when m = 0. Let us define a new operator So, the large- N 'RG vector field': 



+ b{x) 







+ 2<r(x 



9 







_ _ + , where /3 m 

am a;) ora 



m. 



(51) 



<9x^ db(x) 

Then <5oro(-B, E) = 0. The RG vector field may also receive corrections: 5 = 5o + ^ + j^ + -- -. 



4 Small Oscillations around constant classical solutions 

4.1 Constant extrema of N = oo effective action 

Field configurations that extremize the N = oo effective action ro(-B,E) (I41|) dominate the 
path integral over b and a in the saddle point approximation in 1/N. Extrema of To must 
satisfy the large- N 'classical' equations of motion (also called gap equations elsewhere) 

(-V 2 + E)B = and 

_ m 2 } + — JLn(A))s + 0(S - E ) 2 = 0. (52) 

Unlike, say, the gap equations of the non-linear sigma model in the large- N limit, these equations 
are finite and do not require any renormalization, since the divergences in tr log[— V 2 + E] have 



5B 
5E 
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been canceled by the choice of Wo(S) . The complication here is that the equation of motion for 
E has been obtained in a series in powers of £ — S G for constant E D . Let us begin by looking 
for constant extrema (B = B , S = S G ) of Tq. The classical equations of motion become 

Z B = and B 2 Q - m 2 + 2AE = 0. (53) 

Assuming A and m are non-zero, which is the generic situation, there are two types of extrema: 
(S) B a = 0, S G = m 2 /2A where 0(N + 1) symmetry is unbroken and (B) B a = ±m and £„ = 
where 0(N + 1) symmetry is spontaneously broken to 0(N). The vev of the real scalar field 
b must be real. This justifies our choice for the sign of the m 2 a term in Wq(g) in (|40p . In 
the broken phase, the vev Yl — . This is fine for constant backgrounds, since the effective 
potential is entire. But we cannot yet analyze the vicinity of the broken vacuum for non-constant 
backgrounds, since we found the effective action only in an expansion in powers of £ — S D which 
involves inverse powers of E„. Perhaps a cleverer method can be invented to study the phase 
where 0(N + 1) symmetry is broken. There may also be phases where fields are not constant. 
Now, we study oscillations around the vacuum where 0(N + 1) symmetry is unbroken. 



4.2 Mass of long wavelength small oscillations in unbroken phase 



To study oscillations around an extremum of the large- N effective potential, we expand the large- 
./V effective action (|4T1) to quadratic order around a constant background B = B +(3,T, = E -K: 



(54) 



T (B o + /3, E Q + = i J d 4 x (V(3) 2 + |2£ 5 /3 + (B 2 - m 2 + 2A£ ) ? j 



+ < S /3 + 2B [3q + d A 



16vr 2 ' 



n(A) k +■■■ 



We have omitted an additive constant in Tq. Assuming that (B ,T, ) is a constant solution of 
the classical equations of motion (|53p . the linear terms drop out and we get 



r (5 o + ftE o + ? ) 



1 



d 4 x 



{Vf3f 



16vr 2 



tfl(A)? + S ^ 2 + 2B fc + A? 2 + 



(55) 



Here A = — V 2 /S and n(A) is given in (|117p . The lightest 'particles' of the theory are 
the longest wavelength oscillations around the classical vacua. To find the mass of the longest 
wavelength oscillations we may assume that the deviation ? = S — S D is slowly varying in space. 
It suffices to keep the leading term in 11(A) — — -g t ~ W 



A , A 2 A 3 , rp, 
-7T + 9T - "Sn H • Thus 



r (B o + P,E o + <;) 



1 



d 4 x 



(VP) 2 + 



n 



-?(-V 2 )? + S /3 2 + 2B sl3 + A? 2 + 



(56) 



96tt 2 E 

This should describe small amplitude oscillations since we also ignored higher powers of <j. 
In symmetric phase (S) where 0(N + 1) is unbroken, the classical minimum is at B a = 0, 
E Q = m 2 /2A. To expanded around this minimum is 



ro(/?,0 



- / d A x 



^ + 48^^ 2 + ^ 2 + ^ 



(57) 



Recall that a field <fr with Lagrangian (V0) 2 + m 2 (j) 2 has as its longest wavelength excitation, a 
particle of mass m . We deduce that the lightest particle-like excitation of the b field has a mass 
Mb = -j= and transforms in the fundamental representation of 0(N + 1) . Small oscillations of 

a correspond to a particle of mass M a = 4 ^^ m ■ This particle is an 0(N + 1) singlet. There 
could also be heavier particles. 
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5 Leading 1/N corrections due to fluctuations of b and a 



At N = oo, quantum fluctuations of b and a could be ignored, but averaging over them at 
0(1/N) will lead to divergences and violations of scale- invariance in the effective action T . The 
technical difficulty in obtaining r(5,S) = r (B,E) + ^.r x (B,E) is that r (-B,E) ([IT]), when 
expanded, involves arbitrarily high powers and derivatives of E, 



.(58) 



S (b, a) = T (b, (?) = -[ d A x \(Vb) 2 + ab 2 - m 2 a + Xa 2 1 -^\aU(A)a + Oia - a Q 

2 J [ 16tt z I 

II(— V 2 /<t ) is defined in (|117p . The action of our model at this order is 

5(6, a) = S (b,*) + ± J d A xWr (a) + 0(±). (59) 

We must pick Wi(E) so that Ti(B,T,) is finite and scale-free for m = 0. However, this is 
most likely not possible since Wi(E) is independent of B and can at best cancel divergences in 
ri(0, E). The remaining divergences must be canceled via mass and coupling constant renor- 
malization and anomalous dimensions. These will modify the generator of RG flow 

such that r still satisfies the RGE 5F(B, E) = at each order. In order to have a fixed point, 
which would make it natural to have small masses, we will impose (3^ = for at least one value 
of A when m = 0. Since = (150p . the RGE for the effective action at order 1/N reads 

<5 r! + 5^0 = 0. (61) 

5o = —x^-^jt + b-7^ + 2a + rn-J^~ and To ([39]) are known, so we must determine 5\ and Ti . 

5.1 Calculation of effective action at 0(1/N) 

Using (|20f) and doing the gaussian integral at order 1/N we get the change in effective action 

' ■ ■" ' • ' ' • (62) 



tr logrg(B,S) + J d A x Wi(S) 



Tq(B, E) is the hessian of Tq acting on the 2-component column vector (6 a) 



K(B,E)=l V b + S ntJL) x ) (63) 



T± is independent of m even if Tq includes a mass deformation. 



5.1.1 Wi(E) for slowly varying quantum fluctuations 8z general background E(a;) 



To find VFi(E), we calculate T\ for B = and choose Wi(E) to make Ti(0, E) finite and 
scale-free for arbitrary E(x). (However, even for B = 0, E acquires an anomalous dimension, 
see sec. 15.1.21 ) For B = 0, the hessian ([63]) is diagonal, so 



ri(o,E) 



tr log[-V 2 + E] + tr log 



n(-v 2 /E ) 

16vr 2 



+ A 



(64) 
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Since n(— V 2 /E ) is complicated (|117p . we make the further approximation that quantum fluc- 



tuations are slowly varying I. Then we may ignore higher derivatives in n(— V /E G ) and get 



ri(o,E) 



tr log[-V 2 + E] + tr log[-g^- + A] + J Wi(E)d 4 ;c 



(65) 



The 1 st term is identical to what appeared in Tq (|22p and was calculated in an expansion around 
a constant background E G in (l37j) . The 2 nd term can be calculated exactly in zeta function 
regularization. Let ((s) = tr [—V 2 + 96tt 2 XT, ]~ s , then using (f3Tj) we get, 

tr !og[-7^S^ + A] = -C(0)log[96^ 2 S o ]-C'(0) 



where £(0) = 
tr log 



96vr 2 E 
(96vr 2 ) 2 A 2 S 2 



32vr 2 
V 2 



and C / (0) = ^S^log[e" 3 / 2 96. 2 AS o ] 



+ A] 



32vr 2 
(967r 2 ) 2 A 2 S 2 fi 



log[e- 3 / 2 (96^)^AS 2 ]. 



2\2 ■ 



(66) 



96vr 2 E J 32vr 2 
AS cannot be a negative real number. So A = is singular within our approximations. Thus 



ri(o,s) = i 



tr log[-V 2 + E] 



(96tt 2 ) 2 A 2 E 2 Q 

327T 2 



log[e- 3 / 2 (96vr 2 ) 2 AS 2 /M 4 ] + J Wi(E)d 4 x 



.(67) 



The 1 st and 2 nd terms violate scale-invariance and involve a dimensional parameter M as in 
section [3j So Wi(E) is determined by the condition that it must cancel these scale anomalies. 
Aside from the 2 nd term, we had the same condition for Wo(E) m ®. sd3 as in (|40|) . 



W^) = -^{^log 

(96^ 2 ) 2 A 2 E 2 (96tt 2 ) 2 AS 2 . 
2 Si e 3 / 2 M 4 ■ 



2 '-^l + ^ Jfe] ^ - Eo) + ^ log[^](E - E o) 2 



(68) 



Wq + JjW\ is the potential for which To + jjTi is finite and scale-free for backgrounds S(x) 
and B = 0, after including slowly varying quantum fluctuations at 0{\/N). For this choice, 



ri(o,E) 



d 4 x 



16vr 2 



(E - S )n(-V 2 /E )(E - E ) + C>(E - E D ) J 



(69) 



5.1.2 RGE for slowly varying fluctuations and constant backgrounds 

In section [5.1.11 we found Wi(E(x)) that ensures T\ is finite and scale-free for B = 0. But 
there could be further divergences which lead to running couplings and anomalous dimensions. 
We determine (3 m , 7 " and 7 fc , while enforcing /? A = for constant backgrounds B a &: E G . We 
do this by deriving an RGE for the effective action. We assume quantum fluctuations of a are 
slowly varying on the scale of the constant background E Q , so that 



10 Quantum fluctuations need not be slowly varying. We hope to relax this assumption in future work. 

11 We could add a finite term AiE 2 to W\ without violating scale invariance. We did this for Wo, but do not 
for the other W„ , just as we do not add an arbitrary finite counter term X„(f) 4 at n-loop order of A<^> 4 theory. In 
a sense, we allow the most general Wo consistent with the symmetries but the remaining ones are the minimal 
choices that ensure cancelation of divergences and preservation of scale invariance. 
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Here A = -V 2 + E c , D = -V 2 /E D + A and E D = 96vr 2 E . Since B 1 commutes with D, 
detpTg] = detL4L> — B 2 1], i.e., the characteristic polynomial of AD. To calculate it, let 



C(s)= tr [AD- Bl\ 



n 



p 3 dp 



p 4 /E c + (E /E + A)p 2 + (AE D - Bl) 



r- (71) 



The integral converges for $ls > 1 and defines meromorphic £(s) . Changing variables to q = p 2 
and defining 2c = E + AE G , d = (AE — B 2 )T, , £(s) is expressed as 



CM 



q dq 



FIT,* 



16vr 2 Jo (q 2 + 2cq + d) s 32vr 2 



d 



1-s 



c 2 ~ 2s ^(s-i^s + k-d/c 2 



s - 1 



(72) 



£(s) is analytic at s = and we define tr log[AL> — B 2 ] = — C'(0). Then from f)62|) 

ri(s ,s ) = -|c'(o) + ifiWi(s ). 

C'(0) is complicated, but to understand the RGE 5qTi + SiTq = (JBTj) . we only need 



(73) 



C(o) 



32?r- 



(2c 2 -d) 



n 



64?r 2 



E^ + A 2 S 2 + 2B& 



9. 



64vr 2 



{1 + (96tt 2 ) 2 A 2 }E 2 + 192vr 2 5 2 S 



.(74) 



Recall from Sec. E31 that <5 C'(0) = -2C(0). Therefore 
8 Ti = 5 



I C '( ) + inWi(E )) = C(0) + J2Wi(E ). 



(75) 



Thus, the RGE (|6ip for the effective action becomes 



C(o) + -<y nwi(E ) = -<5i r . 



(76) 



From ([35]) we know that Tq(B , E g ) = ^[E D 1? 2 — m 2 E + AS 2 ] . Let us parameterize the leading 
1/N correction to the renormalization group vector field as 



d b d 
1 ^ dm ^ ° dB Q 



7 ff E 



d 
9E C ' 



(77) 



/3 m (A,m) is the beta function 12 ! of m while j b ' a (X,m) are anomalous dimensions. We are 
explicitly imposing the condition /3 A = 0, i.e. that A remains unrenormalized and we have a 
line of fixed points. We only need one fixed point for naturalness, but as the sequel shows, we 
find a line of them within our approximations. Wi(E ), /3 and 7 must satisfy 



C(0) + -<5 JWi(E o ) 



-Sonw l{ E ) + 64yr2 
2 



|(E OJ B 2 -m 2 E + AE 2 



{1 + (96tt 2 ) 2 A 2 }E 2 + 192tt 2 5 2 E 
(-2m(3 m - m 2 7 ' J )E + 2A 7 CT E 2 + (2 7 b + - l a )^ B i 



(78) 



2 P and 7 may receive corrections at 0(1/N 2 ) as part of 82, so f3 m is short for /3] 71 etc. 



19 



We already determined Wi(E) in sec. 15.1.11 from it we calculate 



1 , N ft£ 2 (96vr 2 ) 2 A 2 ft£2 

WWr.ffi,)— (79) 



Putting this in (|78l) we get 

720^ 2 A 2 - 7^)s 2 + 3£ 2 £ = ~ 



64vr 2 ; ° ° 2 



(_ 2m ^ m - mV)£ + 2A 7 CT £ 2 + (2/ + 7 CT )£ D i^ 



• (80) 



This equation must hold for all constant backgrounds B and E . If our model is not renor- 
malizable or does not have a line of fixed points (for slowly varying quantum fluctuations at 
0(1/N)), there would be no choice of /3 m ,j b,cr for which it is identically satisfied. Comparing 
coefficients of monomials in the fields, we get the relations 

2m(3 m + mV = 0, A 7 CT = 720tt 2 A 2 - — ^r, and 7 fc + - 7 a = 3, (81) 

64tt z 2 

whose unique solution i^l (notice that f3 m am so that m = is preserved by RG flow) 

""-Kl^- 360 ^ = 720.= A - y = 3 +I ^-360A (82) 

The resulting RG vector field is 

l + if360vr 2 A yim-^- + { 1 + i- f360vr 2 A - 3 ^rH^o 



A^V 128vr 2 A/J 9m I N \ 1287r 2 A/J °0S o 

+ { 2 4(64k- 720 ' 2A )} E "sl: + O(1/Ar2) - (83> 

For constant backgrounds and slowly varying quantum fluctuations of cr, we have a consistent 
solution of the RGE for which A remains unrenormalized (/3 A = 0). Under these assumptions, 
our model is renormalizable at order 1/N. It remains to see whether renormalizability and 
/3 X = can be maintained (for at least one A), for non-constant backgrounds and rapidly 
varying quantum fluctuations, perhaps via an expansion in inverse powers of S Q . 
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A Examples of naturalness 

By a naturalness explanation for an 'unreasonably small' quantity we mean that the model 
acquires an additional symmetry when that quantity vanishes pQ. In the absence of such a 

13 As noted earlier, A = is a singular limit within our approximations, this is reflected in /3 m ,7 <T ,7 & . 
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symmetry, if the quantity is dimensionless, the 'reasonable' or natural value for it is of order 1 , 
and if it has dimensions, its natural value is of the order of the Planck scale, though there may be 
another scale depending on the context. The symmetry used in a naturalness explanation may 
be either continuous or discrete. Sometimes, it is more convenient to refer to the conservation 
law that follows from the symmetry. The actual small value of the quantity (if non-zero) is 
usually related to explicit breaking of the symmetry and can often be treated perturbatively. 
We give a few examples of naturalness explanations from different branches of physics. It appears 
that this concept explains several small parameters both in experimentally tested theories and 
mathematical models. Indeed, except for the Higgs mass, there is a naturalness explanation 
for every small parameter in the standard model. This gives us some confidence to turn things 
around: if there is an unreasonably small parameter in nature or in a model, then there must 
be some symmetry, which if exact, would make that parameter vanish. In this way, naturalness 
can be a useful guide to model building. 

(1) Small electron mass compared to Planck mass: If m e = QED would gain chiral symmetry. 
The same applies to muon and tau masses. There is a different chiral symmetry for each. This 
puts no constraint on the ratios of lepton masses. 

(2) Small current quark masses: For Nfi avors > 1 if current quark masses are set to zero, QCD 
gains a partial unbroken chiral symmetry SU(Nf)v (isospin for Nf = 2). 

(3) Small coupling constants can be explained by the separate conservation laws for particles, 
gained by setting coupling constants to zero. 

(4) Some near-degeneracies of energy levels in atomic spectroscopy can be explained as due to 
the presence of a symmetry. In hydrogen-like atoms, the difference in energy between levels 
with the same principal and angular momentum quantum numbers n, I but different magnetic 
quantum number m vanishes if we have spherical symmetry. The small energy difference can 
be due to a magnetic field whose direction explicitly breaks spherical symmetry. 

(5) Naturalness in classical mechanics: The fact that planetary orbits are nearly closed and 
nearly lie on a plane, are related to the conservation of angular momentum and the Laplace- 
Runge-Lenz vector in the Kepler problem. Quantum mechanically, in hydrogen-like atoms, the 
'accidental degeneracy' of energy levels with the same value of I is due to a hidden SO{4) 
symmetry whose conserved quantities are angular momentum and Laplace-Runge-Lenz vectors. 

(6) Some near-degeneracies in atomic energy levels can be explained as due to parity symmetry. 
The small splittings are due to parity violation in the weak interactions. 

(7) Experimentally, the mass of a photon is less than 10 -16 eV outside a superconductor. This 
is explained by the exact U(l) gauge symmetry if the photon is massless. 

(8) The near degeneracy in the proton-neutron masses (and pion masses) may be explained as a 
consequence of isospin symmetry. If u and d quarks were degenerate in mass, isospin would be 
an exact symmetry of the strong interactions and the neutron and proton would be degenerate 
in mass (as would the three pions). Isospin is explicitly broken by the quark mass difference as 
well as electromagnetic interactions, which explain the small n—p and ir + — tt° mass differences. 

(9) Pions are naturally light compared to say, p mesons due to chiral symmetry. They are the 
pseudo-goldstone bosons of spontaneously broken chiral symmetry. If the quarks were massless, 
chiral symmetry would be exact at the level of the lagrangian, and be spontaneously broken 
to SU(Nf) v , and the pions would be massless goldstone bosons. But in fact, quarks are not 
massless, this explicit breaking of chiral symmetry gives the pions a small mass. 
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(10) Small neutrino masses: Chiral symmetry for each flavor is exact if neutrinos are massless. 



(11) Parity is an exact symmetry of QCD in the absence of the topological 9 term, which is 
parity odd. Thus, a small QCD theta angle is natural within the theory of strong interactions. 

(12) The scalar field mass in a supersymmetric extension of the standard model can be naturally 
small since if it were zero, the theory would have unbroken global supersymmetry (when the 
super partner fermion is also massless, which would be natural due to chiral symmetry). 

(13) The effective mass (inverse of correlation length) is very small in the neighborhood of a 2 nd 
order phase transition. A naturalness explanation is that at such a transition, when the effective 
mass vanishes, the system gains a new symmetry, scale-invariance. However, one might argue 
that a naturalness explanation is not needed here, since the temperature must be fine-tuned in 
order to have such a phase transition. 

(14) Some linear combinations of correlations in large- iV multi- matrix models vanish because 
of the presence of hidden non-anomalous symmetries [21] . 

(15) It has been suggested [22] that a discrete symmetry that relates real-valued space-time 
coordinates to pure-imaginary ones could ensure a naturally small cosmological constant. 

(16) The amplitudes for gluon scattering either with the same helicity or with only one gluon 
with a different helicity from the others, vanishes in classical Yang-Mills theory. A naturalness 
explanation for this has been suggested, using an effective tree-level supersymmetry |23j . This 
is an example of naturalness in classical field theory. 

(17) GIM mechanism and suppression of flavor-changing neutral currents (FCNC) [17]: The 
small quantity here is the As = 1 strangeness changing neutral current compared to the As = 
neutral current. The GIM mechanism eliminates As = 1 FCNC at tree level by introducing a 
new quark doublet participating in the weak interactions, consisting of charm and the Cabibbo 
rotated strange quark. This is not a naturalness explanation, but could be turned into one by 
specifying a symmetry which ensures absence of FCNC. 

B The quantum effective action 

To be self-contained, we collect a few facts about the effective action (see eg. [H]). The 
generating series for (possibly disconnected) correlations of a scalar field <p is 



The generating series of connected correlations is W(J) = —hlogZ(J). The effective action 
r(<I>) is the generating series of proper vertices (1-particle irreducible diagrams). It is the 
Legendre transform of W( J) : r(<3?) = extj[W(J) — J ■ . The solution of the classical theory 
(tree diagrams) defined by T is equivalent to a solution of the quantum theory of the original 
action S . An implicit integral representation for r(«3>) is 




where 




(84) 




(85) 
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Here $ is the background field and 4> is the fluctuating field. <3? need not solve the classical 
equations of motion, i.e. S"($) need not be zero. To obtain this integral representation we start 
from the relation W(J) = ext$[T(<&) + J ■ <&]. The extremum occurs 

at J = _!EX*). With J 

and $ related this way, the effective action is 

r($) = W(J) - J • $ = -RlogZ(J) - J - $. (86) 
In other words, e~R r ^*^ = Z(J)eR J '*. Inserting the path integral for Z(J), 

e"K r (*) = | [D x ] e -*PW+J-x-J-*]. (87 ) 

Changing variables of integration to = x — $ and using J = J( < I ) ) = — - r ^ , we get the 
implicit formula (|85j) for the effective action. By expanding S(<1> + 0) in powers of (f) and doing 
the gaussian integrals, we get an asymptotic series for in powers of fi (the same formula 

holds with h replaced by 1/iV) 

r($) = 5($) + | logdet S"($) + 0(?i 2 ). (88) 



B.l Large- iV effective potential via dimensional regularization 

We wish to calculate tr log[— V 2 + E ] = /[<i 4 p] log[p 2 + E ] appearing in the large- N effective 
action (j22|) via dimensional regularization. Analytically continuing to n Euclidean dimensions 
and differentiating, we get a convergent integral if n < 2 



Tf vnoTirlinff m a I c 

(2tt)" 



Here [cPp] = ?| f n . Expanding in a Laurent series around n = 4 dimensions, 



t- 5 fe+ E(T - i+ M y w * i) +°t-'>- 

Now ?l tr log[— V 2 + E Q ] = fifLT n . So integrating with respect to E c , 

r 9 , ^E 2 M}E 2 3 , M^E 2 logE „ . . . 

ft tr log[-V 2 + E o] = + ^#(7 - 5 - lo g 4vr) + ^ + + 0(n - 4). (91) 

Integration constant c is independent of E„ and plays no role since it only contributes an additive 
constant to the effective potential. 7 = .577 is Euler's constant. We have a pole part, finite part 
and terms that vanish as n — > 4. Notice that the finite part that transforms inhomogeneously 

under scale transformations lo s s ° j s the same as in cutoff or zeta-function regularization. 

327H 

The choice of Wq that makes the effective action finite and scale-free in the limit n — > 4 is 

ftE 2 fiE 2 logE o , , 

W ° (S °- "> = -16^4) " + AE »' < 92 > 

The finite and scale-free effective action for constant backgrounds is thus 

3/2 -log 4vr)' 

l^ 2 



Fq(B , E g ) = — 



A + — 00.2 ) + X 2 B 



(93) 



We get a line of fixed points parameterized by A , whose definition is scheme dependent 



\ _ A ft( 7 -3/2-log47r) 

^zeta— fn — ^dim— reg t 327T 2 ' 
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C Expansion of tr log[— V 2 + a] in powers and derivatives of a 



C.l Zeta function in terms of the heat kernel 

Let A = — V 2 + a(x) and (a(s) = tr A~ s . Then tr log A = ('(0)- We get an integral 
representation for Ca(s) by making a change of variable 1 1— ► At in the formula for T(s): 

A- s = —— / dt e~ tA If- 1 . (95) 

r(s) Jo 

Now define the evolution operator /i^ = e~ tA which satisfies a generalized heat equation 

^ = = (V 2 -a)h u lim ^ = 1 (96) 

dt t^o+ 

It is convenient to work with the heat kernel ht ip(x) = J dy ht(x,y) ip(y) which satisfies 

dh t (x,y) 



V 2 - a(x) 



h t (x,y) and lim h t (x, y) = S(x - y). (97) 



dt 

Then Qa{s) is the Mellin transform of the trace of the heat kernel: 

] poo I poo p 

U(s) = tr A~ s = —— / dt t 8 - 1 tr e~ tA = —— / dt t 3 ' 1 / d 4 x h t (x,x). (98) 
T(s) Jo r(s) Jo J 

To find h t (x,x) we need to solve (|9"7|) . We find h t (x,y) and take x — > y in the end. For <r = 0, 
ht(x,y) satisfies the diffusion equation and the solution is 

h°(x,y) = J ±^e- L ^ )1 . (99) 



For constant complex a = a Q ([97]) is a PDE with constant coefficients whose solution is 

h ° {x > y) = jhf^ ^^- (100) 

C.2 Short time expansion for heat kernel 

For non-constant a we get an expansion for the heat kernel in derivatives and powers of a for 
small t. We assume that the non-analytic part of the heat kernel is already captured by the 
exact solution (|10U|) . We let a = a a + <j(x) and make the ansatz 

°° e -a t e -(x-y) 2 /4t co 

ht(x, y) = h°(x, y) ^ a n (x, y)t n = — ^ ^ a ™( x ' v)*"- ( 101 ) 

n=0 ' n=0 

The average value of ? need not vanish. However, we assume that V?(x) — > as |x| — > oo so 
that the average value of derivatives of ? and its powers vanish 

J dx (V 2 )V(x) = 0, p,q = 1,2,3... (102) 
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For (jlOip to satisfy the initial condition ([97]) . we need ao = 1. If a is a constant, at = <5o,i- 
In (jlOip . we could absorb e _<Jo * into the infinite series since it is analytic in t, but that would 
amount to throwing away an exact result, so let us not do it. Moreover, e~ CT °* makes the Mellin 
transform (|98p of the heat kernel convergent for $ta > 0, which is necessary to recover Ca(s) . 
The coefficients a n (x,y) are to be obtained by putting (|1U1|) into the heat equation (|9T|) and 
comparing coefficients of common powers of t. We need the expressions 



dt 



dt 



^2 a n t n + h" a n t n 1 and 



o 



i 



V 2 h t = V 2 h° t a n tn + 2V f fc£ J2 a n + h°J2 ^ 



We put these into the generalized heat equation ([571) and get 



dt 



o oo 



^2 a n t n + h° a n t 



n-l 



(103) 



V 2 h° t Y.a n t n + 2V i h° £ t n Vt a n 
o o 

oo oo 

+h° t £ f*V 2 a„ - (a + s)h° t a n t n . (104) 



Using the fact that dth° = (V 2 — a )h° this simplifies to 

oo Y7 LO oo oo oo 

X> + 1) a n+1 1" = 2-^ £ * nv * a « + E rv2 «« - ?E«« tri - 



Now 



« o 



V,; log h° = V,: 



(x - 2/) 2 /4i 



2t 



oo 



So the generalized heat equation becomes 

oo oo oo 

X> + 1) o„+i t n = -(x - y)< X! *" V < a n+i + E PV ' C 


Comparing coefficients of t n determines a n +l given a n with the initial condition ao = 1 
(x - y)i + n + 1 \a n+1 (x, y) = (V 2 - s)a n (x, y). 



Now only a n (x,x) appear in C( s )> so we specialize to 



a n+ i(x,x) 



1 



(n + 1) 



(V 2 - ?)a n (x,x). 



(105) 



(106) 



(107) 



(108) 



(109) 



The first few a n (x,x) are 



ai = -?(x); a 2 = ^(V 2 -0ai = ^(? 2 -V\) 

«3 = J(V 2 -^ = i(V 2 -,)(^-V\) = l(^-^VV-(V 2 f ? ) 
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a 4 = ^(V 2 - ? )a3 = ^(V 2 - ? )(?V\- ? 3 + W-(V 2 )\) 

= ^(v 2 ( ? v 2 ? ) - v 2 , 3 + vV - v 6 , - , 2 V\ + ? 4 - c vV + (110) 

To summarize, the heat kernel expansion is h t (x,x) = ( 6 47rf p Z)o° a nt n ■ If we drop cubic and 
higher powers of then the a n are 



a = 1; ai 



-q a 2 = i(. 2 -V\); «3 = ^V\ + W-(V 2 )\) + 0(a; 
a 4 = ^(V 2 (cV\) + VV- V\ + ? V\) + 0(? 3 ) etc. (Ill) 



For (a{s) we need (a n ) = J d 4 xa n (x,x)/ J d 4 x. Assuming q — ► constant as |x| — > oo and 
V? — > at oo we get (up to terms involving cubic and higher powers of c), 

(oo> = l, <ai> = -(?>, (a 2 ) = ^(? 2 >, (a 3 ) = |( ? V\) + 0(a etc. (112) 

More generally, (a n ) = ^(c(V 2 ) n - 2 c> + 0(c 3 ) for n = 3, 4, 5, ... . 



C.3 Derivative expansion for tr log[— V 2 + a] 

We use ([98]) and the expansion (|10ip to get an expansion for (a(s) in derivatives of a = a Q + q 
1 Z" 00 /" 1 



C(«) 

coo 



r(s) y 
1 



dt / d 4 x t 5 " 1 fct(x,a;) 



E<< 



d£ £ s +™- 3 e - ffo ' 



16vr 2 r(s) ^- jo 

where f2 = J d 4 x. The integral over t is a Gamma function, 

COO 



^a n (x,x)i n 

(113) 







! - - Y(a n )at—T(s + n-2) = a X>Y MEii±_™ " 2 ) 

i67r 2 r(s) v 167r n a o r( 



2-s 



167T 2 



(ao) 



I ( a l) , <02> g(Q3) + l)(a4) 



( s _l)( s _ 2 ) 1 ( s -l)a al al 
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Differentiating and setting s = we get 
C(0) a 2 log a. 







16vr 2 



(«o) (ai) (a 2 ) 



16vr 2 



(ao) 



(ai) 



+ 



16vr 2 



3(a ) (ai) ^ ( n - 3)!(a n ) 



+ E 



log^e" 3 / 2 ] + W(l - loga ) + ^log CTo - Y 



n=3 ° 

oo 



(n - 3)!(a n ) 



n=3 



(115) 



Inserting expressions for (a n ) from (11121) . we get a formula for tr log[— V 2 + a] = —('(0). 

a 2 M 



tr log[-V 2 + a] 



16vr 2 

OO 

-E 

n=3 



1 i r °° i , ( ? ) i r^i , ( ?2 ) i 
2 l0g[ ^ ] + ^ l0g[ T ]+ 2^ l0g<T ° 

(( _ V 2n-4 ?) 



n(ra-l)(n-2)<7? 
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a log[^] , log a a 



16vr 2 



32tt 2 



1 



16vr 2 J 3 n(n - l)(n - 2)(> )™- 2 



+ 0(? 3 ) 



The sum over powers of V 2 can be performed. Let A = — ^- and 

(-A) n A(3A + 2) -2(A + l) 2 log(l + A) 



n(A) = £ 



4A 2 



, n(n + l)(n + 2) 

n=l v ' v ' 

n(A) is analytic at A = 0, 11(A) = -f + |J - ^ + • • •. For large0 A. 

n(A) - 



1, . 3 log A _ . . o. 
--logA + --^- + 0(A- 2 ). 



The final result is 
tr l og [-V 2 + a] = §^log 



=3/2 



+ 



d 1 



X 



16vr 2 



o"o log 



Go 

e 



(116) 



(117) 



(118) 



(119) 



Here cr(x) = a Q + ?(x) , is the volume of space-time and n(A) is defined above. This formula 
is valid for small deviations of a from a constant background a D . The term proportional to ? 
vanishes if a Q is the average value of a. a need not be slowly varying. We assumed that ? 
approaches a constant as |x| — > oo and that V? — > as |x| — > oo. 



Remark: If a is slowly varying, we ignore terms with more than two derivatives to get 



tr log[-V 2 + .] = glog[^] + / 
where a = a a + q and cj is a constant. 



16vr 2 



<r log cr 2 ^V 2 ^ . ~( 3^ 

<r log — k + — ^— ? --7 — + o(s , 

e 2 b(J 



(120) 



C.4 Scale anomaly £(0) for general backgrounds 

Though we only got an asymptotic series for C'(0) around a constant background, we can get 
an exact closed-form expression for its scale anomaly. Under a scale transformation a \— > a 2 a , 



C(s) ' ^ a" 2s C(s) => C'(s) ^ -2C(s)a" 2s log a + a~ 2s C{s) 
C'(0) -> C'(0)-2C(0)loga 

Now we use our result to find the scale anomaly ^(0) : 



121] 



COO K) 



2-s 



n 



16vr 2 



(oo) 



(gi) (02) g(Q3) s(a + l)(a 4 ) 



(«-!)(*- 2) (s-l)(7 a 2 



(122) 



Since all higher order terms are proportional to s, only oq, cl\ and a 2 contribute to C(0) : 



C(o) 



Qa 2 r(ao)_ _ (ai) (02) 



16vr 2 



(To 



erf, 



[i k) k 2 ) 1 

16vr 2 2 o- n 2a 2 



(123) 



Note that —V 2 is a positive operator and er is not a negative real number. 
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D Original potential V(4>) in zeta function regular izat ion 



The simplest way to describe the interactions of our model in the large- N limit is via the finite 
and scale-invariant effective action Tq(B, E) of eqn. (|4ip . Physical intuition and approximation 
methods can be applied to T . The original potential V(<fi 2 /N) is not the effective potential, that 
honor goes to Tq(B,T,) for constant backgrounds. The minima of V{<f> 2 /N) have no physical 
significance. What is more, V((j) 2 /N) is not scale-free, it depends on a scale parameter M . 
But, M is not a parameter of the theory, it is canceled by scale 'anomalies' from quantum 
fluctuations. Moreover, W(o~) is divergent, so V(rj) is not strictly defined independent of a 
regularization scheme. Despite all these warnings, many physicists wish to know what V{rj) is, 
so we find the V(rj) that corresponds to W(a) obtained in zeta function regularization. In the 
large- N limit, we find that for a constant background field a, V(rj) grows as rj 2 / log rj for large 
77 = (j) 2 /N . We haven't yet determined its behavior for small rj. 

Recall that e~ NV ( r, ( x " is the inverse Laplace transform of e - NW i a ( x )) fo r each x: 

^_ -(N/2h)(W(*)+o V ) = ~(N/2h)V( V ) (124) 

We found W(a) only at N = 00, so it makes sense to invert the Laplace transform for large- N . 
We do this here for constant a, for which we have founcf^l W(o~) exactly in ([34 



W(a)+r ) a = r 1 a-m 2 a-^log(Xa) where A = g-^Aft !+3/2]_ (125) 

The branch cut of W(o~) implies C runs from — ioo to ioo avoiding the negative real a axis. 
Putting a = u + iv for u, v £ R, the real and imaginary parts of W(a) + rja = (p + iip are 
(h = l) 

2 \ (n 2 — v 2 ) log AV u 2 + v 2 uv arctan(?j/ii) 

ifi = (ri — m )u k 1 n 

^ w ' 32vr 2 J.67T? 

, . 2 \ (u 2 — v 2 ) arctan(v / u) uv log(X\/u 2 + v 2 ) . . 

^ = ^-rn)v ^ ^ . (126) 

We wee that ?R,(W(a) + rja) — > 00 as v — > ±00 for any u > 0. So the integrand vanishes along 
the lines u±ioo for any u. Thus, the end points of C can be moved from ±ioo to ±ioo + u± for 
any real u± without altering the integral. Simply put, it does not matter along which longitude 
C leaves the south pole or along which longitude it approaches the north pole. 

The general strategy for estimating such integrals is as follows [20J. W(a) + r/a is in general 
complex on C . Its imaginary part ip will lead to a highly oscillatory integral as N — » 00 and 
make it difficult to estimate. The trick is to use analyticity of W(a) + rja to deform the contour 
to a (union of) contour(s) along which Q(W(o~) + rja) is constant and other contours where 
the integrand vanishes. Such contours are called constant phase contours and coincide with the 
steepest contours, those along which the absolute value of the integrand changes fastest. If C is 
such a contour (assumed to be a single one for simplicity), then 

d(T _ e -{N/2h)(W(a)+a V ) = & - ^{W (a)+ar,) j ^ e -(JV/2R)»(W(«r)+<7Tj) _ ( 127 ) 



2m Jc 2iri 



15 We should have a mass scale M to set the scale for the logarithms. We set M = 1 in this section. 
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Now we have eliminated the oscillating phase and for large- N, the asymptotic behavior is 
determined by the local minima of p> = 3?(W(<t) + err])) along C. Since cp diverges at the end 
points of C , local minima must occur at interior points of C . Moreover, there must be an odd 
number 1, 3, 5, • • • of such local minima along C . At any one, the directional derivatives of both 
<p and ip vanish in the direction tangent to the curve. Since (p + iip is analytic, it follows that 
these local minima of ip are saddle points, i.e. d a (W(a) + arf) = 0, where two or more steepest 
curves intersect. Not all saddle points of W(a) + err/ need lie on C and those that don't will not 
contribute to the asymptotic behavior of the integral. 

Suppose a = cr s (v) is the only saddle point along the constant phase contour C. The 
integrand attains a local maximum at a s along C and decays exponentially in either direction 
away from a s . The contour can be approximated by a straight line tangent to C at a s and of 
length e on either side. p>(a) = 3?(W + rja) is approximated by its quadratic Taylor polynomial 
around a s , whose linear term vanishes. Now we let e — > oo. <p(cr s ) gives the leading contribution 
while the quadratic term in its Taylor series gives a gaussian integral proportional to -7= . So 

e -fv{ v ) = e -f ^(<7 3 ) e -f^K)J_ 0( 1 
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V( V ) = (p (a a )+iiP(a a ) + o(^-\ (128) 

If there is more than one saddle point on C , we add up their contributions in this formula for 
V{rj). Moreover, if the saddle point a s is on the real axis, then ip(cr s ) = does not contribute. 

So our job is to find a convenient constant phase contour and identify the saddle points on it. 
In practice, we look for saddle points and then a suitable contour. The saddle point condition 
for W(a) + arj for given 77, m 2 and A is 



16vr 



2 



(r]-m 2 ) = oTog(Aoye). (129) 



Taking imaginary and real parts, it is a pair of transcendental equations (for A = 1 and h = 1) 



— + — log (u +v) = — u arctan(w/u) and 

9 . 9 , U 111 _ . q O v / i \ / \ 

167T (77 — m ) = — + — log (u +v ) — v arctan(u/u). (130) 

The fundamental domain for the arctangent is taken as — ir < arctan < tt. We must solve 
for a = u + iv assuming it is off the negative real axis. Any u > for v = satisfies 
the first condition (i.e. saddle points can lie on the positive real a axis), but there are other 
possibilities. The imaginary part of the saddle point condition is also satisfied along a curve 
(found numerically) in the u—v plane that encircles the origin and is symmetric under reflections 
about either axia^£| and lies within the rectangle^ |u| < e _3//2 , \ v\ < 4j. However, the second 

saddle point condition is satisfied on this curve only for a limited range of values of r\ — m 2 , 
namely m 2 > 7? — m 2 > — no 1 7= for — e -3 / 2 < u < and — 00 1 7= < 77 — m? < —m 2 for 

J c — / — ?,2-K^/e — — 327r v 'e — ' — c 

< u < e~ 3 / 2 where m 2 = he %~ . For 77 — m 2 in this range, saddle points could occur on the 
above mentioned curve as well as the positive real a axis, making their analysis more involved. 



16 This follows from the even and oddness of the condition in u and v respectively. 

17 The limiting values are obtained by solving the first saddle point condition for small v and u respectively. 
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For now, we set aside the behavior of V(rj) for small 18 ! rj, i.e. rj — m 2 < m 2 . For rj — m 2 > m 2 
the only possible saddle points are located on the positive real a axis (v = 0) . In this case, the 
positions of the saddle points are given by the solutions to equation (|129f) where a = u is real. 
Since we assumed r\ — m 2 > m 2 , the LHS is in particular positive, and there is a unique solution 
o s which can be found recursively 

V fj fj 

a s = = = = = = = • • • (131) 

log(AVeo-s) i g ( hfm \ log(A-v/e?7) - log log(X^/ea s ) 



\ log(A y/ea s ) 

where fj = 167r 2 (r/ — m 2 )/Ti. Thus, for sufficiently large rj — m 2 , there is only a single saddle 
point a s . Moreover, we have checked numerically that there is a constant phase (actually zero 
phase, since tp(a s ) = 0) contour that starts at the south pole, passes through a s and goes to 
the north pole. For large rj — m? , the leading approximation for the position of the saddle point 
is (with A = 16-7T 2 yfeX/U) 

16-k 2 (rj — m 2 ) o 

as ~* ¥1 — rTT 2vi as V-™ ^oo. (132) 

n,log[A(?7 — m z )J 

The original potential is given by V(rj) = W(a s (rj)) + rja s (rj) + C(log N/N) . Using the saddle 
point equation (1129p . we simplify this to 

V{rj) = W(a s ) +rja s = ^ + ~a s (r, - m 2 ). (133) 

For large rj — m 2 we geJ^I 

8ir 2 (rj — m 2 ) 2 



V{r,) 



h log[A(?7 — m 2 )] 



21og[A(r] — m 2 )] 



as rj — m — > oo (134) 



Recalling that rj = we see that up to a multiplicative factor, V((f) 2 /N) ~ ^^y^ for large 

cj) 2 /N and fixed m. Thus, the original potential grows logarithmically slower that a quartic 
potential in the large- N limit. It would be interesting to find the behavior of V(rj) for small rj. 
This requires a careful study of saddle points and constant phase contours for small rj . 
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